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Recently we studied the B, Bj, D and D^ meson decay constants using various treatments for the 
heavy quark. For B mesons, we determined fe, fss^ and fsj fa with NRQCD bottom quarks. We 
then combined the ratio fg^ / fg and another very precise determination from HPQCD for fs, using 
heavy HISQ quarks, and extracted /b with 2% total errors. We also calculated fo, fo,, and fojfo 
using HISQ charm quarks. Here we review our results and briefly discuss their implications for 
the determination of the CKM matrix elements \Vcci\ and \Vcs\- 
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1. Introduction 

Investigating the flavor structure of the Standard Model (SM) is important for its own sake; 
however, it is even more interesting since it can lead to physics beyond the SM. Furthermore, the 
data accumulations and new analysis emerging from the LHC suggest more and more that the 
Higgs is very close to the SM Higgs. They even do not find any hint of new particles yet. In this 
situation, precise understanding of the SM in the flavor sector becomes more critical. 

Decay constants of heavy-light pseudoscalar mesons have been studied from lattice QCD for 
quite some time. We can determine the corresponding CKM matrix elements by combining decay 
constants from theory and decay rates from experiments. Moreover, decay constants are basic 
quantities related to many hadronic quantities. For instance, /g is an important input parameter for 
inclusive determinations of the CKM matrix elements. The decay constants can also be used for 
testing the lattice formalism, since the calculations have typically smaller errors and the procedures 
are relatively straightforward. 

This paper presents new calculations for B, B^, D and Dg meson decay constants from HPQCD. 
We have completed the projects, and published the results in two papers [|l]][|2|]. So, essentially this 
proceeding consists of a brief summary of the two papers and a discussion of their impact on the 
CKM matrix elements \Vcd \ and \Vcs\- 

2. B and Bs meson decay constants 

We used the MILC AsqTad Nf = 2 + 1 gauge configurations with NRQCD (Nonr^elativistic 
QCD) b quarks for this project. The previous HPQCD calculation [^] used AsqTad light and 
strange valence quarks; however, in this work we used the highly improved staggered quark (HISQ) 
action for the valence quarks. The HISQ action has much smaller discretization effects, so one can 
expect improvements in the continuum extrapolation errors. We include one more ensemble (FO); 
40^ X 96 with mi/m, = 0.0031/0.031, which is a more-chiral fine ensemble. Details of the lattice 
configurations are in Tab. |l]. 

Moreover, this calculation is done with the new scale parameter ri values. HPQCD was using 
ri = 0.321 (5) fm extracted from T splittings [Q]. In 2010 HPQCD published a much more accurate 
ri determination, ri = 0.3133(23) fm, based on several physical quantities and an improved con- 
tinuum extrapolation with 5 lattice spacings [^. We needed to re-tune valence quark masses due to 
the scale changes, and updating for such different new settings is one of the main purposes of this 
analysis. We used the spin averaged T mass to tune the bare bottom quai^k mass, and (fictitious) rjs 
for the strange quark mass. Fig. |I| shows the tuning for b (left) and s (right) quarks. As one can 
see, our physical target meson mass has a large eiTor compared to the deviations of the alignment 
of the tuning measurements. We found that it is important to tune quark masses precisely up to 
the statistical errors or ri/t? eiTors of the tuning measurements. This precise tuning ensures correct 
estimation for of chiral and continuum extrapolation. If one has large deviations between the 
tuning measurements more than its statistical eiTors, then the chiral and continuum extrapolations 
may suffer from additional systematic eiTors. 

We calculated operator matching factors in full QCD at one-loop through order a,, 
aas, and a,^^^. These matching calculations were presented separately at this conference [^. We 
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Set 


ri/a 


m/ /nis (sea) 


Nconf 


Ntsrc 


L^xNt 


CI 


2.647 


0.005/0.050 


1200 


2 


24^ X 64 


C2 


2.618 


0.010/0.050 


1200 


2 


20^ X 64 


C3 


2.644 


0.020/0.050 


600 


2 


20^ X 64 


FO 


3.695 


0.0031/0.031 


600 


4 


40^ X 96 


Fl 


3.699 


0.0062/0.031 


1200 


4 


28^ X 96 


F2 


3.712 


0.0124/0.031 


600 


4 


28^ X 96 



Table 1: Simulation details on three "coarse" and three "fine" MILC AsqTad ensembles. Ncnf is the 
number of the configurations that were used in the simulation, and Ntsrc is the number of time sources per 
configuration. 
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Figure 1: Tuning of the b (left) and s (right) quark masses. 



also employed random- wall sources for the HISQ propagators and Gaussian smearing sources for 
the NRQCD propagators. 

Including all statistical and systematic errors, we obtained 

/B=0.191(9)GeV, /b, =0.228(10)GeV, (2.1) 

and 

^ = 1.188(18). (2.2) 

JB 

These calculations are adefinite improvement on our previous calculations fs = 0.190(13)GeV, 
fs^ = 0.231(15)GeV, and fsj fs = 1.226(26). The largest source of enms is the operator match- 
ing error for the decay constants. The ratio fsjls has very small errors, since most of the matching 
factors are canceled. 

If one calculates the decay constants without matching factors, one could reduce around 5 % 
errors down to 1 ~ 2 % errors. Recently, HPQCD has calculated fs^ without matching factors |^], 
and obtained fs, = 0.225(4)GeV with only 1.8 % en^ors. Essentially, this very precise calculation 
utilizes the HISQ action for the b quark. The HISQ action can be used to simulate the charm quark. 
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Figure 2: Comparisons of results for fg (left) and fB, (right) from this analysis, previous HPQCD, Fermi- 
lab/MILC and ETM collaborations. 

but for the bottom quark it is very difficult with current technology. The clever idea was that in 
fact we can simulate a heavy quaik heavier than the charm quark but lighter than the bottom quai^k. 
Once one gets heavy meson coiTelators depending on multiple heavy quark masses, then one can 
extrapolate to the physical bottom quark mass. In this way, one can determine the decay constants 
without matching factors. 

Of course, we can apply this heavy HISQ method for fs, but it would be difficult. First of 
all, the light quark is much more expensive than the strange quark, and for we need to perform 
additional chiral extrapolation. Thus, for now we fix fs by combining fs^ / fB from the NRQCD 
analysis and fs, from the heavy HISQ analysis; 

-1 

X ^ o.l89(4)GeV. (2.3) 

NRQCD 

This fs result with 2 % total error is the most accurate fs available today. Comparisons of results 
for fs (left) and fs, (right) are shown in Fig. ^ 

3. D and Ds meson decay constants 

With the same simulation setting shown in Tab. [T[ we determined D and D., meson decay con- 
stants. The discretization error of the HISQ action starts at i^{as{amhYv^ /c^) and ^((am/,)^v^/c^), 
and this provides enough accuracy to simulate relativistic charm quarks on current typical lattices. 
So, we apply the HISQ action for all valence quarks including the charm quark. Thus, we can 
evaluate the decay constants without matching factors, since the HISQ action exhibits the chiral 
symmetry in the continuum limit. The decay constant can be written with the heavy-light axial 
vector current = ^qYnY5^c, with q = s or: d, d& 

<0\A^\D>=p^fD,, (3.1) 
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Oldri 


New ri 




241(3) 


246(4), 248(3) 


fn 


207(4) 


208(3) 


foJfD 


1.164(11) 


1.187(12) 




231(15) 


228(10), 225(4) 


fB 


190(13) 


191(9), 189(4) 


IbJIb 


1.226(26) 


1.188(18) 


Ik 


157(2) 


159(2) 


fn 


132(2) 


132(2) 



Table 2: Decay constants from HPQCD with the old ri = 0.321 (5)fm and the new ri = 0.3133(23)fm. The 
unit of the decay constants is MeV, and the ratios are dimensionless. 



We can express the decay constant in terms of the pseudoscalar density PS = ^qjs^c, as we used 
for light meson decay constants, and /k, 

fDq = ^TIT^ < O^^P? > • (3-2) 

As we did for the B decay constants, we also re-tuned the charm quark mass for the new 
r\ = 0.3133(23) fm. We used T]^ mass to fix the charm quark mass. HPQCD updated fj)^ with the 
new ri in 2010 [|8|] already, so this work is mainly to update fn with the new r\. Our final results 
are 

/z) = 208.3(1. 0),f«a3.3),v..MeV, /o, = 246.0(0.7),,,,. (3-5).vs.MeV, (3.3) 

and, 

^ = 1.187(4),,„J12),,,., (3.4) 

which show good agreement with our previous determinations. 

One interesting question would be what the impact of the scale change is. Tab. ^ summarizes 
HPQCD 's determinations of the decay constants with the old and new scale factor r\. We found 
no significant effect due to the scale change. (One exception is for fr,^ with the old r\ and the most 
accurate result with the new r\) It appears that for the decay constants re-tuning of quark masses 
largely compensates the shift from overall scale change. Thus, predicting the impact of the scale 
change before the actual calculations would be risky. We will investigate the impact of the scale 
change further in the future. This study would lead to better estimation of systematic errors for the 
scale setting. 



In Fig. 0, we compare fo and fo, results from FNAL/MILC [|9p, HPQCD [§|][110|], ETMC \^ 
and PACS-CS [||]. The comparisons include FNAL/MILC's prehminary results with Nf = 2 + 
1 + 1 including simulations at the physical pion mass, and ETMC's preliminary results with Nf = 
2+1 + 1. In their preliminary results, they achieve a good precision that is comparable to our best 
results, and they show very good agreement with HPQCD. 

Combining our decay constant results with branching fractions from experiments, we can ob- 
tain corresponding CKM matrix elements. See Fig. ^ for the comparisons. For \Vcd\, one can 
immediately notice that the leptonic determination and the semileptonic determination are in good 
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Figure 3: Comparisons of results for fo (left) and fo^ (right). The results of this proceeding are shown 
under HPQCD 2012. 
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Figure 4: Comparisons of results for \Vcd\ (left) and \Vcs\ (right) from leptonic decays and semileptonic 
decays. 



agreement. Those two results were obtained with completely different systematics for both lattice 
and experiment analysis, since the decay channels are quite different. This is a highly non-trivial 
check for lattice formulation and experiments. Those two lattice determinations of \Vcd\ demon- 
strate good agreement with the unitarity point as well. Thus, we do not see any signature of new 
physics here yet. We note that now the precision of the lattice determination of \Vcd\, especially 
with the preliminary branching fraction result from BES III is actually better than the accuracy 
of the determination from neutrino experiments. So far, the PDG quotes the neutrino experiment 
result for \Vcd\- This is simply because lattice determinations had much larger eiTors in the past. 

For I V„ I , the situation is more interesting. The right plot of Fig. ^ shows that the semileptonic 
determination and the unitarity point are in good agreement. However, for the leptonic determina- 
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tions, it shows some discrepancies depending on the decay channels and experiments. If we only 
consider average leptonic determinations with averages of HFAG, then the |Vc.v| represents more 
than 1 a discrepancy between the unitarity point, which is the identical to observation of the fo. 



puzzle. The /^^ puzzle [ |13[ ] was a 4 a tension between experiments and lattice determinations of 
in around 2010. The puzzle is no longer a puzzle, since new experiments and lattice analysis 
results have moved closer to each other. Now the difference is about 1.6 a. When experiments 
determine fo,, they use the unitarity |- Thus, the f^, puzzle indicates a difference between iVc^l 
from the unitarity point and our determination from the leptonic decay of Ds meson. In Fig. ^ we 
show the results with two Belle's 2012 preliminary results, the first is for Ds — )• pLV and the second 
is for Ds TV, these results suggest that we need to wait to see the experiments attain more accu- 
racy. It may still be possible to see some discrepancy in determinations of |Vc.v|. and this could be a 
hint for new physics. 
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